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Abstract We report on the analysis of XMM-A''ewto;j observations of three G-type stars in very different evolutionary phases: 
the weak-lined T Tauri star HD 283572, the Zero Age Main Sequence star EK Dra and the Hertzsprung-gap giant star 31 
Com. They all have high X-ray luminosity (~ 10^' ergs"' for HD 283572 and 31 Com and ~ 10^^" ergs"' for EK Dra). We 
compare the Emission Measure Distributions (EMDs) of these active coronal sources, derived from high-resolution XMM- 
Newton grating spectra, as well as the pattern of elemental abundances vs. First Ionization Potential (FIP). We also perform 
time-resolved spectroscopy of a flare detected by XMM from EK Dra. We interpret the observed EMDs as the result of the 
emission of ensembles of magnetically confined loop-like structures with different apex temperatures. Our analysis indicates 
that the coronae of HD 283572 and 31 Com are very similar in terms of dominant coronal magnetic structures, in spite of 
differences in the evolutionary phase, surface gravity and metallicity. In the case of EK Dra the distribution appears to be 
slightly flatter than in the previous two cases, although the peak temperature is similar. 
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1. Introduction 

During the last decade, the analysis of high-resolution X-ray 
spectra of late-type stars, o btained with EUVE, XMM- 
Newton and Chandra (e .g. iMpnsignori Fossi et all Il995t 
.Schmitt et al. 1996: Gu del et alJ Il997t iGriffiths & Jordari 
19981: iLaming & Drak3 Il999t ISanz-Forcada et alJ l2002t 
Argiroffi etalJl2003l) . revealed that the thermal structure of 
coronal plasmas is better described by a continous Emission 
Measure Distribution, EMD, rather than by the combination 
of a few isothermal components usually employed to fit low- 
and medium-resolution spectra. Since the coronal plasma is 
optically thin, the EMD of the whole stellar corona can be 
viewed as the sum of the emission measure distributions of 
all the loop-like structures where the plasma is magnetically 
confined; therefore, it can be used to derive information about 
the pr operties of the co ronal structures and the loop popula- 
tions jPeres et alj|200ll) . In particular, the studies mentioned 
above have indicated that the coronae of intermediate and 
high activity stars appear to be more isothermal than coronae 
of solar-type stars, and that the bulk of the plasma emission 
measure is around log T ~ 6.6 for stars of intermediate 
activity and up to log T ~ 1 for very active stars. The latter 
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result is consistent with the one previously obtained from the 
analyses of Einstein and ROSAT data, i.e. that there is a good 
correlation between the effective coronal tempe rature and the 
X-ray emission level (see, for example, Schm itt et alJll990l 
|PreibisclJl9'97l) . 

The observation that in active stars a considerable amount 
of plasma steadily resides at very high temperatures, which are 
achieved on the Sun only during flaring events, led to the hy- 
pothesis that a superposition of unresolved flares may heat the 
plasma causing an enhanc ed quasi-quie scent coronal emission 
level. Following this idea. iGiidell ll 1 997h showed that the time- 
averaged EMD resulting from hydrodynamic simulations of a 
statistical set of flares, distributed in total energy as a power 
law, could be made quite similar to the EMD of stars of dif- 
ferent activity level (and age). In particular, he obtained dis- 
tributions with two peaks and a minimum around lOMK; the 
amount of the hottest plasma (at ~ 12-30 MK) decreases with 
decreasing Lx (or, equivalently, with increasing age) and, at the 
same time, the first peak moves towards lower temperatures. 

A qualitatively different scenario for the evolution of the 
EMD with acti vity has been p roposed by J. Drake (see Fig. 2 
in the review by iBowver et al. 2000): the distribution increases 
monotonically from the minimum, which occurs at log T be- 
tween ~ 5 and 6, up to the peak at coronal temperatures; the 
location of the peak shifts towards higher and higher tempera- 



tures (up to log 7 ~ 7 in the most active stars) ror increasing 
X-ray activity level. Along with the shift of the peak, the steep- 
ness of the ascending part of the distribution increases. 

In the pictures sketched above, the shape of the EMD 
changes with the stellar activity level; however, it is not yet 
understood which stellar parameters (luminosity, surface flux, 
surface gravity, evolutionary phase, or others) have a major 
role in determining the physical characteristics of the domi- 
nant coronal structures and, hence, the properties of the whole 
Emission Measure Distribution. 

In order to investigate this issue, we have examined the 
cases of three G-type stars, in different evolutionary phases: 
the Pre-Main Sequence star HD 283572, the Zero-Age Main 
Sequence star EK Draconis (HD 129333) and the Hertzsprung- 
gap giant star 31 Com (HD 111812). Here we report on the 
analyses of recent XMM-Newton observations of these bright 
targets, characterized by similar and relatively high X-ray lumi- 
nosities (Lx ~ lO-'^'ergs"' for EK Dra, and Lx ~ 10^^ ergs"' 
for HD 28357 2 and 31 Com) wit h respect to the Su n. Previous 
analyses (e.g. IC.iidel et al. 1997:.Avres et al.. 1998- Favata eHll 
119981) showed that the characteristic coronal temperatures of 
the stars of our sample lie around 10^ K; the EPIC and RGS 
detectors on board XMM are very sensitive to this temperature 
regime, allowing us to get rather accurate and reliable informa- 
tion about the plasma Emission Measure Distributions of these 
stars. 

The a nalysis of th e XMM observation of 3 1 Com was re- 
ported in IScelsi etakl ll2004l) . while the reconstruction of the 
EMD of HD 283572 using a high-resolution spectrum is pre- 
sented here for the first time. For ease of comparison with these 
two sources, we have also re-analyzed the XMM observation of 
EK Dra using the same method employed for 3 1 Com and HD 
283572, thus ensuring homogeneity of the results; note how- 
ever that independent analyses of the same XMM obser vation 
of EK Dra have been pub lished since 2002 (e.g. Gude l et alJ 
l2002HTelleschi et al.l2003l) and more recently and comprehen- 
sively byEHssEh^^LllE^) in the context of a study of so- 
lar analogs at dififerent ages. The latter work is complementary 
to our present study because it considers stars having similar 
mass, size and gravity, but largely diff'erent Lx and coronal tem- 
perature. 

This paper is organized as follows: we describe the three 
targets in Sect. Hand we present the observations in Sect.|3] In 
Sect.0]we describe the data reduction and the methods used for 
the analyses of EPIC and RGS spectra. The results are shown 
in Sect. Island discussed in Sect.|6l 

2. The sample 

In Fig.[2we plot the positions of the sample stars in the H-R 
diagram, to show their respective evolutionary phases. We used 
visual magnitudes, B-V color indexes and distances measured 
by Hipparcos; we assumed negligible optical extinction in the 
cases of EK Dra and 3 1 Com, coherent with the low interstellar 
absorption used in the analysis of their X-ra y spectra (Sect . |5j, 
while we used a visual extinction Ay = 0.57 istrom et alJl989l) 
and Eb-v ~ Av/3 for HD 283572. 
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Figure 1. Positions of HD 283572, EK Dra and 31 Com in 
the H-R diagram. We have superimposed pre-main-sequence 
(dashed lines) and post-main-sequence (solid lines) tracks for 
the mass va lues reported in the plot . The evolutionary models 
are those of lVenturaetaD ( Il998albl) . except for th e 2Mp) pre- 
M.S. track, for which we have used the model of ISiess et alJ 
(2000). All tracks are calculated for solar photospheric abun- 
dances. 



The latter star is a member of the Taurus-Auriga star 
forming region and its age is estimated to be ~ 2 x lO^yr 
(IWalter etalJll988l) . HD 283572 shows no sign of accretion 
from a circumstellar disk, which characterizes the earlier stage 
of classical T Tauri stars; the decoupling from the disk al- 
lowed this star to spin up cosiderably, due to its contraction, 
up to several tens of kms"' (vsin/ = 78kms"', see Table 0, 
probably with a consequently enhanced dynamo action and a 
very high X-ray luminosity (Lx ~ 10^' ergs"'). From Fig.^ 
we deduce that HD 283572 will be an A-type star during its 
main sequence phase, and we estimate a mass between ~ 1 .5 
an d ~ 2Mo, in agreement wit h the estimate of 1.8 + 0.2 Mq 
bv lStrassmeier & Rice* ('l998a'). The radius of HD 283572 has 
been derived by Walter et al. ( 1987) through the Barnes-Evans 
relation, R ~ 3.3Ro at an assumed distance of 160pc, which 
becomes 2.7 R© at the new distance of 128 pc measured by 
Hipparcos; more recently, Strassmeier & Rice (1998a) com- 
bined photometric measurements, rotational broadening and 
Doppler imaging technique to determine the radius of HD 
283572 in the range 3.1 - 4.7 Rq, with a best value of 4.1 Rq. 
Due to the uncertainties of these estimates, we decided to con- 
sider both of them. We anticipate that our main results are only 
weakly affected by the choice of one of these values. 

EK Dra is a G1.5-type star with mass and radius about 
equal to the solar values. It has just arrived on the main se- 



lable 1. Stellar parameters. Distances are measured by Hipparcos; (0.3 - 8 keV) are derived in this work rrom i - 1 models; 
gravities are determined from the corresponding M and R, and surface fluxes from the corresponding Lx and R. In the last column, 
the references for Lboi are indicated in the entries. 
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quence, thus representing an analog of the young Sun. Beca use 
of its age (~ 7 x 10^ yr, LSoderblom & Clem ents 1987), it 
suffered littl e magnetic brak ing and its short rotational period 
(~ 2.7 days, iGuinan et alJ20 03) makes it a bright X-ray source 
(Lx~ lO^^ergs-i). 

The more massive (M ~ 3 Mo) giant star 3 1 Com (age 
~ 4 X 10** yr, Friel & Boesgaard 1992) has already evolved out 
of the main sequence and now it is crossing the Hertzsprung- 
gap. The position of 3 1 Com in the H-R diagram and the evo- 
lutionary models indicate a spectral type late-B/early-A on the 
main sequence; therefore, this star has developed a convective 
subphotospheric layer and a dynamo only in its current post- 
main sequence evolutionary phase (Pizzolato et al. 2000). The 
X-ray luminosity is ~ 7 x lO-'" erg s"' . 

The stellar parameters of the three targets, with the relevant 
references, are summarized in Tabled For HD 283572 we re- 
port both estimates, mentioned above, of the stellar radius and 
the corresponding values of gravity and surface X-ray flux. 

HD 283572, EK Dra and 31 Com were chosen because 
their stellar parameters offer the possibility to get useful in- 
sight into their coronal properties from the comparison of their 
EMD. Note, in particular, that while the X-ray luminosity of 
31 Com and HD 283572 are about equal and larger than that 
of EK Dra by about an order of magnitude, EK Dra and HD 
283572 are the stars with the highest surface fluxes, whose val- 
ues exceed significantly that of 3 1 Com, by about one order of 
magnitude. Note also that the different evolutionary phases im- 
ply different stellar internal structures; moreover, these targets 
have quite different gravities, implying different pressure scale 
heights and possible changes in the properties of the dominant 
coronal loops. 

Finally, the rapidly rotating stars HD 283572 and 31 Com 
are putative single sources: this avoids difficulties in the inter- 
pretation of the results, due both to the uncertain origin of the 
emission, in case of multiple components, and to the possibility 
of an enhanced activity as found, for example, in tidally-locked 
RS CVn systems. On the contrary, EK Dra has a distant com- 



panion jDuquennov et ^ljl99ll). whose mass is likely between 
0.37 and 0.45 jCiidel et al Jl 995ah . lCTiidel et alJjl 995W) 

found that the X-ray and radio emissions are modulated with 
the rotational period, strongly suggesting that the coronal emis- 
sion comes predominantly from the G star If we assume that 
the secondary star has M ~ 0.4 and age ~ 70 Myr, and has 
a saturated corona (the worst case), its X-ray luminosity would 
be ~ 10^^ erg s"', so we might expect contamination of the X- 
ray emission of the G star from the companion at most at ~ 10 
% level. 

3. Observations 

The observations of HD 283572, EK Dra and 3 1 Com were per- 
formed with XMM-Newton respectively on September, 5, 2000 
(PI: R. Pallavicini), on December, 30, 2000 (PI: A. Brinkman) 
and on January 9, 2001 (PI: Ph. Gondoin). The non-dispersive 
CCD cameras (EP IC MOS and EPIC pn. Turner et al. lOOlli 
IStriider et al. Il200lh . lying in the focal plane of the X-ray tele- 
scopes, have spectral resolution 7? - E/AE ~ 5-50 in the range 
0.1 - 10 keV, while the two reflection grating spectrometers 
(RGS. Iden Herder et alJ200lh provide resolution R ~ 70-500 
in the wavelength range 5 - 38 A (0.32 - 2.5 keV). 

For the present study, we considered only the EPIC pn and 
RGS data; in Table |2l we report details on the instrument con- 
figurations and on the observations. 

At the time of these observations, both CCD 7 of RGS I and 
CCD 4 of RGS2 were not operating. These CCDs correspond to 
the spectral regions containing the He-like triplets of neon and 
oxygen, respectively. Note also that the RGS 1 spectrum of HD 
283572 is entirely missing, due to instrument setup problems in 
the early phase of XMM-Newton observations; hence we have 
no information on the O vii triplet for this source. 

4. Data analysis 

We used SAS version 5.3.3, together with the calibration files 
available at the time of the analysis (June 2002), to reduce the 



lable 2. Log or the XmM-Newton observations. 
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" Exposure time for the analysis of the quiescent emission (Q.E.), i.e. excluding the time intervals affected by proton flares, 
occurred in the cases of HD 283572 and 31 Com, and by the source flare in the case of EK Dra (see Sect.lO. 

* Mean count-rate in the 1.2 - 62 A (0.2 - lOkeV) band for pn and in the 5 - 38 A (0.32 - 2.5 keV) band for RGS (1st order 
spectrum) relevant to the Q.E. Exposure. 



data of HD 283572 and 31 Com; the data of EK Dra were re- 
duced with SAS version 5.4 and the analysis was performed in 
September 2003. We generated all pn responses with the SAS 
RMFGEN and ARFGEN tasks. 

Good Time Intervals were selected by excluding those time 
intervals showing the presence of presumable proton flares in 
the background light curve extracted from CCD 9 of the RGS, 
following den Herder (,2002): we cut the intervals where the 
count-rate exceeds 0.1 ctss"' for 31 Com and 1.6 ctss"' in 
the case of HD 283572 (whose observation is contaminated by 
high level of background), while we did not exclude any inter- 
val in the case of EK Dra. 

In order to obtain X-ray light curves and spectra, we ex- 
tracted the events from a circular region (~ 50" radius) within 
CCD 4 for HD 283572, while we used annular regions (~ 
7.5" - 50" radii) for 31 Com and EK Dra, because the relevant 
data were affected by pile-up. In all cases, background photons 
were extracted from the rest of CCD 4, excluding the sources 
and their out-of-time events. 

4.1. Light curves 

Figure |2l shows the pn background-subtracted light curves of 
the sources, with a 200 s time binning. The light curve of 31 
Com is the only one that is consist ent with the hypothesis of a 
constant emission (see Sect. 3.1 in IScelsi et'aill2004i) . 

In the case of HD 283572, there is evidence of variability 
of the emission, on a time-scale of the order of 30 ks, which 
is not a typical flare event. The reduced is 5.9 (229 d.o.f.) 
in the null hypothesis of a constant emission; the variability 
amplitude, calculated as 0.5 [max(rate)-min(rate)]/mean(rate), 
is ~ 20%. There is a less pronounced variability on a time-scale 
of ~ 10 ks. From Tables^and^ we note that the duration of the 
observation is about one third of the stellar rotational period, 
hence a large fraction of the stellar surface was visible during 
the pointing; this suggests that at least part of the variability is 
due to an inhomogeneous distributiuon of active regions over 
the stellar surface. 

The light curve of EK Dra clearly shows the presence of a 
flare; the vertical lines in the figure mark the start and the end of 
the flare, obtained as the minimum and maximum times where 
the hardness-ratio,///? = (H - S)/{H + SY , systematically ex- 

' We have evaluated the soft emission count-rate, 5, in the 0.3 - 
1 keV band and the hard emission count-rate, H, in the 1 - lOkeV 
band. 



ceeds by more than 1 cr the average HR value calculated from 
time intervals before and after the flare. We excluded the time 
interval of the flare from the emission measure analysis (Sect. 
[4.31 . since we want to study the thermal properties of the qui- 
escent corona, and we analyzed the flare separately. The qui- 
escent emission of EK Dra is still variable, yielding a reduced 
- 6.2 (190 d.o.f.) against the null hypothesis of a constant 
source; the variability is on a time-scale of ~ 15 ks and its am- 
plitude (calculated as above) is ~ 16%. 

4.2. EPIC PN spectra 

We have performed global fitting of the EPIC pn spectra (Fig.|3} 
with the aim of deriving, from multi-component thermal mod- 
els, the initial guess of the continuum level for the line measure- 
ments in the RGS spectra. Moreover, the abundances of some 
elements (Si, S, Ar, Ca) can be better determined from pn spec- 
tra, rather than from RGS spectra, thanks to the wider spectral 
range of the former - which includes the strong K-shell lines of 
the relevant H-like and He-like ions - and to the higher photon 
counting statistics-. 

We analyzed these spectra with XSPEC vll.2 and we 
found that an absorbed, optically-thin plasma with three 
isothermal components provides an acceptable description 
of each of them (see results in Sect. 15. U . The models 
are based on the Astrophysical Plasma Emission Database 
(APED/ATOMDB VI. 2) and have variable abundances; we 
adopted the criterion of leaving free to vary only the abun- 
dances of those elements (O, Ne, Mg, Si, S, Fe, Ni, in some 
cases Ca and Ar) with strong and clearly detectable line com- 
plexes in EPIC spectra. The abundances of the other elements 
were tied to that of iron, their best-fit values being poorly con- 
strained when left free to vary. 

We eventually used the high-energy tail of the pn spectrum 
also to check the high-temperature tail of the emission measure 
distributions, as described in the next section and in Appendix 

m 

Finally, we performed time-resolved spectroscopy of the pn 
data of EK Dra during the flare, to get information on the prop- 

- The Sixiii-xiv lines fall also in the RGS spectral range, but the 
statistics are usually very low and the calibration of the effective area is 
less precise at these wavelengths; nonetheless, the results of our anal- 
ysis show that the Si abundances derived from pn data are consistent 
with those obtained from RGS spectra within statistical uncertainties 
(see Table|4j. 
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Figures. EPIC pn spectra of HD 283572, EK Dra and 31 Com with their best-fit model spectra (the parameters of the models 
are listed in Table|3. The spectra of 31 Com and HD 283572, with their relevant best-fit models, have been shifted by -0.5 and 
-1-0.5 dex for clarity. 



erties (in p articular the size) of the flaring loop, employing the 
method bv lReale'e t al. ( 1997). This analysis and its results are 
reported in AppendixiBl 

4.3. Emission Measure Reconstruction 

The approach we adopted for the line-based an alvsis of the 
RGS s pectra of each star is discussed in detail in IScelsi etaP 
( I2OO 4') together with a study of its accuracy; here we limit our- 
selves to report the main points of our iterative method. 

We employed the software package PlNTofALE 
jKashvap & Drake 2000) and, in part, also XSPEC, and 
used the APED /ATOMDB VI. 2 database which includes the 
[Mazzott a et al.l {^98) ionization equilibrium. 

We first rebinned and co-added the background-subtracted 
RGS 1 and RGS2 spectra for the identification of the strongest 
emission lines and the measurement of their fluxes. In this 
latter step, we adopted a Lorentzian line profile and we as- 
sumed initially the continuum level evaluated from the 3-T 
model best fitting the pn spectrum, because the wide line wings 
make it impossible to determine the true source continuum be- 
low ~ 17 A directly from the RGS data, in particular in the 
~ 10 - 17 A range, where the spectrum is dominated by many 
strong overlapping lines. Then, with the aim to reconstruct the 
Emission Measure Distribution (EMD) vs. Temperature, we 



selected a set of lines, among the identified ones, with reliable 
flux measurements and theoretical emissivities. Most of them 
are blended with other lines, so the measured spectral feature 
is actually the sum of the contributions of a number of atomic 
transitions; accordingly, we evaluated the "effective emissiv- 
ity" of each line blend as the sum of the emissivities of the lines 
which mostly contribute to that spectral feature. Moreover, we 
carefully selected only iron lines not blended with lines of other 
elements, because the procedure we employed (see below) uses 
these iron lines in the first step of the EMD analysis, and esti- 
mates of the abundances of the other elements are not yet avail- 
able at this step. 

We performed the EMD reconstruction with the Markov- 
Chain Monte Cai-lo (MCMC) method by Kashyap & Drakq 
(1199811) . This method yields a volume emission measure dis- 
tribution, EM(Tk) - dem{T It) AlogT, and related statistical 
uncertainties AEM{Tk), where demij) - WgdV/dlogr is the 
differential emission measure of an optically thin plasma and 
A log T = 0.1 is a constant bin size; the method also provides 
estimates of element abundances, relative to iron, with their sta- 
tistical uncertainties. The iron abundance is estimated by scal- 
ing the emission measure distribution assuming different metal- 
licities and by comparing the synthetic spectrum with the ob- 
served one at /I > 20 A in the RGS spectrum (this is a spectral 
region free of strong overlapping emission lines). Finally, we 
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are important tests ror the reliability or the amount or plasma 
in the high-temperature tail of the EMD. 

We also checked the consistency between the continuum 
level assumed for flux measurements and the continuum pre- 
dicted by the EMD. In fact, since our method is iterative, the 
continuum assumed for flux measurements in the RGS range 
is adjusted at each iteration for consistency with the EMD, 
and it may become different from that predicted by the 3-7" 
model best-fitting the pn spectrum, which is adopted as initial 
guess. Therefore, this procedure ensures that possible cross- 
calibration offsets between pn and RGS do not aff'ect the final 
EMD. 




31 Com 
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Figure 2. Background-subtracted pn light curves of HD 283572 
(upper), EK Dra (middle) and 31 Com (lower), in the 0.2 - 
lOkeV band and with time bins of 200s. The vertical lines in 
the light curve of EK Dra mark the time interval of the flare, 
excluded from the analysis of the quiescent emission. 



checked the solution obtained with the MCMC by comparing 
(i) the line fluxes predicted from our solution with the measured 
ones and (ii) the pn model spectrum, based on the reconstructed 
EMD, with the observed pn spectrum at £ > 2keV. These 
checks are illustrated respectively in Fig. |6l and in Appendix 
IaI taking the case of EK Dra as an example (similar results 
were obtained for the other two stars). In particular, the correct 
prediction of the O vii-viii line fluxes allowed us to check the 
reliability of the amount of plasma in the low-temperature tail 
of the EMD; analogously, the correct prediction of the Fe xxm- 
XXIV line fluxes and of the high-energy tail of the pn spectrum 



5. Results 



5.1. 3-T models 



Figure |3l shows the pn spectra with their relevant 3-T mod- 
els, obtained by fitting the data in the 0.3 - 8keV range. The 
best-fit parameters of the models are listed in Table|3l 

The presence of the Fexxv 6.7 keV emission line in all 
these spectra is indicative of hot coronae, as expected from ear- 
lier works and confirmed by our analysis. Note the large best-fit 
EM values of the hottest components for all stars, comparable 
to the EMs of the cooler components; in particular, the hottest 
plasma dominates the corona of HD 283572, as also con firmed 
by the analysis of Chandra spectra jAudard et alJ2004 . 

The line complexes of Mgxi-xii (~ 1.3 - 1.5 keV), Sixiii- 
XIV (~ 1.8 - 2.1 keV) and S xv (~ 2.5 keV), as well as the large 
bump between 0.6 and 1 keV due to the Fe xvii-xxiv, Ni xix-xx 
and Ne ix-x lines allowed us to constrain the abundances of Mg, 
Si, S, Fe, Ne and Ni. These complexes are less evident in the 
spectrum of HD 283572, as a consequence of the significantly 
lower metallicity with respect to the other two stars; instead, the 
Ca XIX line complex (~ 3.9 keV) is most prominent in the spec- 
trum of this star (Fig.O and the estimated abundance of this el- 
ement is higher than for the other two cases. Note also that we 
were able to constrain the Ar abundance for EK Dra, thanks to 
the clearly visible lines of Ar xvii at ~ 3.1 ke V. In the other two 
cases we linked the abundances of Ca and/or Ar to that of Fe 
assum ing the same ratios as in the solar corona jGrevesse et alJ 
Il992h . Moreover, we used the results of the EMD analyses to 
fix the coronal C/Fe abundance ratio for EK Dra and 3 1 Com to 
respectively 0.7 and 0.25 solar, and the N/Fe ratio for EK Dra 
to 0.5 solar 

Finally, we could not constrain the interstellar absorption in 
the directions of EK Dra and 3 1 Com with the fitting procedure, 
so we fixed them at the relatively low valu es of 3 x 10 ' ^ cm~^ 
and 10^^ cm~^ mea s ured, respectively, by iGiidel et alJ ( 1 19971) 
and IPiskunov et alJ ( Il997 ^. On the contrary, the spectrum of 
HD 283572 is significantly absorbed, as expected from the lo- 
cation of this star in the Taurus-Auriga star forming region. The 
hydrogen column density we derived from the fit is compatible 



lable 3. Best-ht models or the bFlC pn data (in the U.3 - s keV band), with 90% statistical conhdence ranges computed ror one 
intere sting parameter at a t ime; nominal errors on T, and EM, are at the 10% level. Element abundances are relative to the solar 
ones jGrevesse et aljfl993) . Mean temperatures are calculated as < T >= ^Li 2?=i EMj. Abundances and hydrogen 

column densities without errors were fixed as explained in the text. 
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Figure 4. Co-added RGS spectra of HD 283572 (upper), EK Dra (middle) and 3 1 Com (lower) with the identification of the most 
prominent lines; the bin size is 0.02 A. 



with Ay and consistent with previous res ults obtained by fitt ing in all cases, while O vii, Mg xi-xii and Si xm-xiv emission lines 
ASCA, ROSAT, Einstein and SAX data (Favata et al.'l92Sl)- 
5.2. Emission Measure Distributions and abundances 



The rebinned and co-added RGS spectra (Fig. |3 show emis- 
sion lines from Fe xvn-xxiv, Ne ix-x, O viii and Ni xix-xx ions 



are visible only tor bK Ura and 31 Com, and the N vii line in 
the case of EK Dra only. Actually, we could not identify any 
line outside the wavelength range 10 - 20 A in the spectrum of 
HD 283572, because of contamination from high background 
and lack of the RGS 1 spectrum altogether. The reconstruction 
of the EMD of EK Dra and 3 1 Com was based on about 40 
lines, while we used 25 lines in the case of HD 283572 as a 
consequence of the lower quality of its spectrum. 

The derived EMDs are plotted in Fig.|5l note that the algo- 
rithm we used is not able to constrain statistically the values of 
the emission measure in all the temperature bins. We show in 
Fig.|6lthe observed-to-predicted fluxes for the case of EK Dra, 
which is representative of the spread of these ratios obtained 
in our analyses, and in Fig.0we compare the observed spec- 
tra and the model spectra generated with the solutions (EMD 
and abundances) found in this work. The elemental abundances 
are shown in Tabled) we estimated the iron abundances (rela- 
tive to the solar value) of HD 283572, EK Dra and 31 Com at 
0.7 + 0.2, 1.2 + 0.2 and 1.4 + 0.2, respectively. Table |5]reports 
the ratios^ R - f li and G — (f + i)/r relative to the O vii triplet, 
and the estimates of electron temperatures, densities and pres- 
sures, averaged over the region w here the triplet forms, using 
the theoretical curves bv lSmith et al. (2001). 

In Fig. |8] we show the element- to-iron abundance ratios 
for the three stars, ordering the elements for increasing First 
Ionization Potential (FIP). Whenever both EPIC pn- and RGS- 
based estimates were available, we always report the latter in 
the plot, because we consider the values derived with the RGS 
the most accurate. It is worth noting that, despite widely dif- 
fering methods employed to derive elemental abundances, we 
obtained consistency between the pn- and RGS-derived abun- 
dance ratios, except for Ne in the case of 31 Com (the value 
indicated by the pn is two times larger than the RGS one) and 
for Ni in the case of EK Dra and HD 283572 (the values ob- 
tained with the PN are larger than the RGS ones by factors of 2 
and 4 respectively). 

The patterns of abundances vs. FIP are similar in the cases 
of 3 1 Com and HD 283572, with an initial decrease (with re- 
spect to solar photospheric values) down to a minimum around 
carbon, followed by increasing abundances for elements with 
higher FIP (> 1 1 eV). This pattern is also similar to what was 
found for the young active star AB Dor bv lSanz-Forcada et alJ 
( 2003 ), but it is less evident in the case of EK Dra. Note that 3 1 
Com and EK Dra have iron abundances differing from that of 
HD 283572 by about a factor of 2, hence the pattern of abun- 
dances vs. FIP appears to be almost independent of the global 
coronal metallicity. 



6. Discussion 

XMM-Newton data allowed us to derive the plasma emission 
measure distributions for our three targets and their coronal el- 
emental abundances; in particular, the EMD of HD 283572 has 
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Figure 5. Distributions of emission measure derived from RGS 
data. Values without error bars are not statistically constrained 
by the MCMC algorithm. Note the different ordinate scale in 
the plot of EK Dra with respect to those of HD 283572 and 3 1 
Com. 



been derived here for the first time using a high-resolution spec- 
trum. Our results are sufficiently well determined and homoge- 
neous for the purpose of a detailed comparison of the coronal 
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Figure 7. Model spectra compared to the original RGS spectra. 



lable 4. Ratios between elemental and iron coronal abundances, relative to the solar photospheric ratios tCjrevesse et alJll992l) . 
derived from RGS data; errors are at 68% confidence level. For completeness, we also report the absolute iron abundance. The 
PN-derived values are shown for purpose of comparison. For each star, the number of lines used for the EMD reconstruction is 
reported (the number of lines of a given element is shown in parenthesis near the relevant (RGS) abundance value). 
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Figure 6. Comparison between observed fluxes and the fluxes 
predicted with the EMD model, for lines used in the EM re- 
construction of EK Dra; Fe: open diamonds, Ne: triangles, 
Mg: open squares. Si: filled diamond, Ni: filled circle, O: filled 
squares, N: asterisk, C: open circle. 




□ HD 283572- 
A EK DRA 
O 31 COM 



Co Ni Mg SI S C N Ar Ne 

Figure 8. Element-to-ir on abundance ratios, r elative to the so- 
lar photospheric values d Grevesse etaD 'l992V for HD 283572 
(squares), EK Dra (triangles) and 31 Com (diamonds). The el- 
ements are ordered by increasing FIR 



properties of the selected stars. We recall that these stars are 
in different evolutionary stages, but share the characteristic of 
being active (high X-ray luminosity) G-type stars. Our analysis 
has confirmed that the three stars have very hot coronae, with 
similar average temperatures (~ 11 - 12 MK for EK Dra and 
31 Com, and ~ 16 MK for HD 283572). 

A remarkable result of this work is the close similarity of 
the emission measure distributions of HD 283572 and 3 1 Com, 
which have similar Lx as well. Both distributions have a well- 
defined peak at Tp = 10^ K and, in the range log T ~ 6.5 - 7, 



Table 5. Pressure estimates with O vii. 
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they are proportional to ~ T^, where the exponent of the power 
law has a formal confidence interval between ~ 3.4 and ~ 6.6; 
there are also indications of a significant amount of plasma at 
temperatures hotter than Tp (up to log T ~ 7.6) and, at least 
in the case of 31 Com, in the range log T ~ 6 - 6.2. We re- 
call that we are not able to statistically constrain the emission 
measure in all the temperature bins, and hence to get informa- 
tion on the exact shape of the distributions below log T ~ 6.5 
and above log T ~ 7.3, yet the presence in both stars of a non- 
negligible amount of plasma up to log T ~ 1.6 has been veri- 
fied, as described in Sect. 14. 31 through the correct prediction of 
the Fe xxm-xxiv line fluxes and by comparison with the high- 
energy tail of the observed EPIC spectra (App. |A}, while the 
correct prediction of the O vii-vm lines allowed us to verify the 
presence of cool plasma down to log T ~ 6 in the EMD of 
3 1 Com (the O vii line is not available in the spectrum of HD 
283572). Note, also, that the shape of the constrained part of the 
distribution of 31 Com and the presence of significant emission 



measure at log 1 ~ 6 - 6 .2 suggest that a minimum in the EMD 
of this star occurs around log T ~ 6.4-6.5, while some caution 
is needed for the case of HD 283572. 

The EMD of EK Dra is, on average, about one order of 
magnitude lower than the two previous ones. The distribu- 
tion has a maximum at log T - 6.9, with a more gradual de- 
crease towards high er T than in the two previous cases. Using 
ASCA/EUVE data, ICiiidel et alJ (Il997l) derived for EK Dra 
an EMD essentially bimodal, with two significant peaks near 
7 MK and 1 8 MK. While we also find little plasma at tempera- 
tures below ~ 3MK and our value of Tp (~ 8 MK) is roughly 
consistent with their first peak, we do not find either a deep 
minimum around lOMK or strong evidence for a second max- 
imum at 18MK. 

It is not straightforward to derive a low-T slope for this dis- 
tribution, essentially due to the secondary peak at log T - 6.5. 
If we exclude this temperature bin, on the grounds that it be- 
longs to a cooler population of coronal structures (see below), 
a fitting in the range log T - 6.6-6.9, yields a slope of 5.3 + 1.7: 
this is compatible, within errors, with the slopes derived for HD 
283572 and 3 1 Com. Instead, if we consider the secondary peak 
as a fluctuation produced by the emission measure analysis, the 
slope in the range log T - 6.5 - 6.9 turns out to be 3.0 + 1.2. 
To address the issue of the EMD slope, we have investigated if 
a solution smoother than the one presented in Fig.|5]can give a 
good description of the observed line fluxes as well. Assuming 
a distribution with its ascending part (log T < 6.9) proportional 
to and with the high-temperature tail (log T > 6.9) decreas- 
ing as T^^, we searched for the minimum on the subset of 
the iron lines, by varing the exponent jS and a global renor- 
malization factor. The minimum is obtained for - 1.3 
(the relevant EMD is shown in Fig. |9] together with the recon- 
structed one), and the observed fluxes for all the selected lines 
are quite acceptably reproduced (Fig.^| to be compared with 
Fig. |6|l. Note also that we used, for this test, the relative el- 
ement abundances reported in Table 0] However this solution 
over-predicts the Ovm 18.97 A/Ovii 21.60A line ratio: while 
the observed ratio of line counts is /?o,obs = 10.0 + 1.7 (1 cr 
error), the ratio predicted by the smooth solution is Rq ~ 16, 
against Rq - 10.6 predicted by the MCMC solution. This line 
ratio does not depend on the oxygen abundance, but it is es- 
pecially sensitive to the shape of the distribution at low tem- 
peratures (log T < 6.8), hence its predicted value could be 
lowered by a small enhancement of the emission measure at 
logT = 6.0 - 6.3 without affecting significantly the fluxes of 
the iron lines"*. Interestingl y, our smo oth solution is similar to 
that found recently by Telleschi et aL C2004) using the same 
XMM data, but a different inversion method. Thus, we are not 
able to get strong constraints on the low-T slope of the EMD 
of EK Dra and further investigation is needed. 

We want to interpret the shape of the bulk of the emission 
measure distributions, around 10^ K, in terms of loop struc- 
tures. First, we note that, at this high coronal temperature, the 
pressure scale height Hp in each of the three stars is of the or- 
der of the corresponding radius (Table|6|l. We assume that the 
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Figure 10. Comparison between observed fluxes and the fluxes 
predicted with the smoothed EMD (dashed line in Fig. |9}, for 
lines used in the EMD analysis of EK Dra; Fe: open diamonds, 
Ne: triangles, Mg: open squares. Si: filled diamond, Ni: filled 
circle, O: filled squares, N: asterisk, C: open circle. 



structures responsible for the observed emission of each star 
have characteristic sizes smaller than the rel evant pressure scale 
height. In fact, for 31 Com we showed in 'Scels i et alJ ( l2004 
that loop heigths larger than the stellar radius would be hardly 
compatible with the absence of variability in the emission of 
this star, because of the very low filling factor implied by this 
solution; in Appendix IbI we analyze the flare on EK Dra and 
conclude that the height of the flaring loop, which likely be- 
longs to the family of structures dominating the X-ray emis- 
sion, is a few 10'" cm; finally, Favata et al. (2001) analyzed a 
flare on HD 283572 and showed that the size of the involved 
structure is ~ 0.3 Rf Under this hypothesis, the pressure is ap- 



proximately unirorm inside each loop, implying that the emis- 
sion measure distribution of a si ngle loop depends only on 
its maximum temperature r^ax (fMaggio & PereJl99(Tl) . with a 
functional form EM(T) cc T" for T < T^ax, with a = 3/2 in the 
case of loops with constant cross-section and uniform heating. 
Considering that the EMD of the whole stellar corona is the 
sum of the EM{T) of individual loops, the total EMD would be 
proportio nal to T" for T < m in(rmax); hence, following the ap- 
proach bv lPeres et alJ ll200ll) . we interpret the constrained part 
of the EMDs of HD 283572 and 3 1 Com as due to a population 
of loops, each of them having EM{T) cc T^, since the EMDs 
of these stars are approximatively power-laws (with exponent 
~ 5) in the temperature range mentioned above. Consequently, 
the simplest interpretation we derive from the comparison of 
the emission measure distributions of HD 283572 and 3 1 Com 
(Fig- HI see also Fig. below) is that the coronae of these 
stars are very similar in terms of dominant coronal structures, 
in spite of their different evolutionary phases and gravities, as 
well as coronal abundances. Since this latter parameter plays 
an important role in the energy balance through the radiative 
losses, we might expect that different abundances result in dif- 
ferent temperature and density profiles along a loop, and hence 
in different coronal EMDs; however, this is not the case in these 
two stars. Moreover, we stress that HD 283572 and 31 Com 
show similar EMDs in spite of the difference in X-ray surface 
flux by about one order of magnitude. In conclusion, we infer 
that all the above parameters have only a minor role in deter- 
mining the properties of the EMD of these stars, which instead 
appear to be mainly determined by the high and nearly identical 
X-ray luminosity. 

The high index (~ 5) of the power law which best approx- 
imates the ascending part of the EMDs of these stars also sug- 
gests that the dominant coronal loops of very bright sources 
(with Lx ~ 10^' ergs"') may have different properties to the 
solar ones. In such stars, the physical processes that lead to 
emission measure distributions significantly steeper than those 
observed in low-luminosity stars, such as the Sun or a Cen 
tD r a ke et al. 19 97). still remain to be understood; as discussed 
in Iscelsi et alJ 120041) . a possible interpretation of such steep 
slopes, which should also characterize the EM(T) of the sin- 
gle structures (see above), is that the heatin g of the coro- 
nal loops is located mai nly at their footpoints jTesta & PeresI 
l2003t iTesta etaljEooi ). Another possibility are expanding 
loops, although quite extreme expansion factors are needed 
dSchriiver et al J 19891) : also, both effects might be at work. 

The interpretation of the EMD of EK Dra is more uncer- 
tain. If we assume that the distribution reconstructed in this 
work and shown in Fig.|5]is a good approximation of the actual 
EMD of this star, then a family of hot loops with EM(T) oc 
appears to be responsible for the bulk of the distribution around 
10^ K, while the emission measure at logT < 6.6 may be 
due to a cooler family of loops whose properties we are not 
able to investigate with the available data, because of the lim- 
ited information on the low-temperature plasma; otherwise, 
if a smoother distribution applies (EMD oc T^, Fig. |9}, the 
structures dominating this corona might have a less steep pro- 
file of the emission measure vs. temperature with respect to 
the cases of HD 283572 and 31 Com, yet steeper than the 



lable 6. Stellar radii and pressure scale heigths calculated tor 
T ~ 10^ K. 
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T^^^ slope characterizing uniformly-heated loops with constant 
cross-section. 

In Fig. the EMDs of the three studied stars are shown 
together with the emission mea sure distribution of the Sun 
dPeres et alj|200 0') and ^ Bootis ("Laming & Drakd l 19991) . flie 
latter being a G-type star of intermediate activity, with Lx ~ 
lO^*^ erg s ' . While the EMD of the Sun peaks at Tp ~ 10^-^ K, 
with an ascending part (log T ~ 5.7 - 6.2) proportional to T^^^, 
and shows no significant amount of plasma at temperatures 
above ~ 10^'^ K, the EMD of ^ Bootis is intermediate between 
those of the Sun and our active stars, in terms both of Tp and 
of the overall amount of emitting plasma, as well as with re- 
gard to the steepness of the EMD preceding its peak. Actually, 
^ Bootis is a double star (G8-I-K4V), but the X-ray emission 
is thought to be dominated by the primary G8 (Schmitt 19il, 
and we have included it in this picture because it is one of the 
very few stars of intermediate activity whose EMD has been 
reconstructed from high-resolution spectra. 

The comparison shown in Fig. ^2 between the EMDs of 
stars with increasing luminosity, going from the quiet Sun to 
HD 283572, suggests a transition in the steepness of the dis- 
tribution, which, in turn, may reflect changes of the properties 
of the dominant coronal loops. Fig.^Jis in agreement with the 
h ypothesis of increasin g steepness with increasing Lx reported 
in lBowver et al.l ll200dl) : on the other hand, their picture is not 
reflected completely by our results which do not indicate, at 
coronal temperatures, a monotonic increase of the EMD up to 
its peak, at least in the cases of the bright star 31 Com and pos- 
sibly also of HD 283572; in fact, the XMM data available for 
these two stars suggest the presence of plasma at T ~ 10* K 
and a minimum of the EMD around T ~ 10^-^ K. 

Finally, it is important to note that the sample is still limited, 
and new observations are required to have a more complete 
scenario, as well as to study in greater detail the "cool" and 
"hot" tails of the emission measure distributions of very active 
stars. In this respect, the recent work by Telleschi et al. (2004) 
provides complementary results which may help to bridge the 
gap between solar-type stars and stars with very high activity 
levels. 

Appendix A: Checking the high temperature tail of 

the EMDs 

As stated at the end of Sect. 14.31 we checked the reliability 
of the hot tails of the EMDs, not constrained by the MCMC 
method, by comparing the model spectrum with the high- 
energy (E > 2 keV) PN spectrum, which is more sensitive to 
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Figure 11. Emission mea sure d istributions of the coronal plasma of the quiet Sun jPeres et alJ2000t from Yohkoh/SXT data), ^ 
Bootis {Laming & Drak^ l 19991 from ASCA/EUVE data), EK Dra, 31 Com and HD 283572 (this work). For each of the three 
active stars, only the constrained part of the EMD is shown. 



very hot plasma (it contains, in particular, the complex around 
6.7 keV largely dominated by Fexxv and several spectral re- 
gions dominated by the bremsstrahlung continuum). This com- 
parison is shown in Fig. lA.ll for the case of EK Dra, similar 
results have been obtained for HD 283572 and 31 Com. In this 
plot, the solid line is the pn model spectrum derived from the 
reconstructed EMD, whose hottest part is reported (with the 
solid line) in the inset. In the high-energy tail of the spectrum 
(2-8 keV) the model fits well the data (xl ~ 0.98, 88 d.o.f.), 
indicating reliability of the presence of a sizeable amount of 
plasma in the high-temperature tail of the EMD. We have also 
examined both the cases of a flatter hot tail with a total emis- 
sion measure twice as large as the previous case (approxima- 
tively proportional to T^^l^, shown by the dashed line in the 
inset), and of a tail which decreases as T^^ (dotted line), with a 
total emission measure four-fold lower than the first case. The 
corresponding model spectra give a poor fit to the data {xl ~ 4 
and 1.4, respectively, consider that P(xl > 1.4) ~ 0.8%). 

Appendix B: Analysis of the flare on EK Dra 

The flare observed on EK Dra (Fig. IS has a duration of ~ lOks 
and is characterized by a rather wide (~ 3 ks) peak at ~ 2 cts/s, 
subtracting the mean count rate of the quiescent phase, and a 
secondary maximum at ~ 0.5 cts/s which follows the initial de- 
cay phase. For the analysis of this flare, we performed time- 
resolved spectroscopy of the EPIC pn data and employed the 




10"^ I . . . . . 

7.0 7.2 7.4 7.6 



, loq T (K) 
10-* I I . . . . . 

1 

Energy (keV) 

Figure A.l. EPIC pn spectrum of EK Dra (crosses) and the 
model spectrum derived from the reconstructed EMD (solid 
line); the dashed line is the model spectrum obtained in the 
case of a distribution equal to the reconstructed one up to 
log T - 7.3, but with a flatter hot tail drawn with the dashed 
line in the inset; the model spectrum shown with the dotted 
line, instead, is relevant to a steeper decrease of the hot tail. 



approach bv lReale et alJ lll997h to derive the size of the flaring 
loop, assuming that the impulsively heated plasma was con- 



tined in a single structure. We rerer to that paper ror a detailed 
explanation of this method, and to Reale et al. (2004) for its ap- 
plication to a flare (on Proxima Centauri) observed by XMM- 
Newton. 
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Figure B.l. The upper panel shows the light curve of the flare 
on EK Dra, in the 0.3 - 10 keV band of the EPIC pn, obtained 
by subtracting the average quiescent emission from the total 
light curve in Fig.|2 The vertical lines mark the phases of the 
flare. The middle and the lower panels show the time evolu- 
tion of the temperature and the emission measure, respectively. 
Filled symbols are relevant to the model with variable global 
metallicity, open symbols to the model with z - 0.83 solar 



tions with time of the temperature and the emission measure 
are shown in the middle and lower panels of Fig. IB. II 

Unfortunately, the results of our analysis for this specific 
event are aff'ected by several uncertainties, hence the loop size 
we obtain is to be taken with caution. 

The half-length L of the loop is a function of the e-folding 
time of the light curve, Tlc, of the maximum temperature of 
the flaring plasma, Tobs.max, and of the slope ( of the trajectory 
(during the decay) in the density-temperature diagram. To de- 
rive tlc we need the light curve of the flare alone, which is the 
total light curve minus the quiescent emission; yet, in our case, 
the latter component is comparable to the flaring one and it is 
not constant, as akeady shown in Sect. 14.11 We approximated 
the quiescent emission with its constant mean value in order to 
obtain the light curve of the flare, shown in the upper panel of 
Fig. IB. II This curve is not characterized by a well-defined peak 
followed by an exponential decay, therefore it is not possible to 
derive an accurate e-folding time from it. We estimated tlc 
by evaluating the time, after the peak, when the count-rate has 
fallen down to 1 /e times the maximum value and we obtained 
Tlc ~ 4000 s. 

Another difliculty arises from the second maximum^ of the 
flare, which "breaks" the decaying phase. Since the statistic is 
not very high for such a kind of time-resolved spectral analysis, 
we were able to derive the trajectory in the - T diagram 
during the decay and its slope only from two points (segments 
d\ and dl). The slope ^ is estimated to be ~ 1.2 using the 
results of the fittings with variable z, and ~ 1 in the second 
case; these values indicate that sustained heating was present 
during the decay of the flare^. 

The maximum observed temperature was evaluated to be 
T^obs ~ 4.2 X 10^ K in both cases (segment r in Table IB. it : we 
conclude that the size of the flaring loop is of the order of a few 
10"^ cm, i.e. substantially smaller than the stellar radius. 

Finally, we derived the average temperature Tgq of the flar- 
ing plasma at equilibrium, i.e. after the flare has totally decayed 
and the loop returned to its quiescent conditions, by fitting the 
temperature values of segments r, pi, p2 and dl with the func- 
tion r = A e '^^ -I- Teq. Although the errors are rather large, we 
obtained a minimum of for = 6 X 10^ K, in the case of 
variable z, and T^q = 5 x 10^ K, in t he case of fixed z; f rom these 
two values we obtain (see Eq. 4 in iReale et al.l2004l) the max- 
imum temperature of the loop in quiescent conditions (reached 
at its apex): r^ax.eq = 8.4 x 10*'K and r^ax.eq = 6.8 x lO^'K, 
respectively. These estimates are around the peak temperature 
of the EMD of EK Dra and may indicate that the loop where 
the flare occurred belonged to the family of loops contributing 
to the bulk of the X-ray emission of this star. 



We divided the observation during the flare into 6 segments 
(Fig. IB.lt . so as to have ~ 4000 - 5000 counts in the pn spec- 
trum of each of them. These spectra were fitted in XSPEC us- 
ing, as a model, the fixed 3-T model in Table|3] describing the 
quiescent emission, plus a fourth component, which gives the 
temperature and the emission measure of the flaring plasma. 
We made the fittings both using a variable global metallicity 
(z) and fixing it to the quiescent value (z = 0.83 solar). The 
results of the fittings are reported in Table IB. II and the evolu- 
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